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ABSTRACT p 


Observations of the temporal evolution of the wave-formed ripples are analyzed 
in terms of geometric properties, migration rate, and forcing by wave and current 
velocities. Three weeks of bedform observations were obtained using underwater video 
of a sheet of laser light projected on the bed at the Monterey Bay Inner Shelf Observatory 
in thirteen-meter depth water. The bed consists of fine sand (mean grain size 0.12mm). 
Low to moderate narrow-band swell waves occurred during the observation period. 
Ripple geometry consisted of orbital and suborbital vortex wave ripples, and relic ripples 
left after larger wave events. The bedforms generally changed size and shape when the 
grain roughness Shields parameter exceeded a critical value of 0.04 with the exception of 


the last event bringing into question the appropriate value of 0 


,,1n nature. Ripple 
migration was offshore at rates of 2-10 cm/day during active sediment transport events. 
Skewness of velocities (low and high passed) were calculated to explain offshore ripple 
migration, but showed no preferred direction. Mean currents were weak and also showed 
no preferred direction. Significant positive (offshore) correlation was obtained between 
the short-wave envelope and infragravity waves indicating wave-group forced bound long 


waves (surf beat) combined with stirring by the short waves might explain the offshore 


sediment transport and ripple migration. 
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1. INTRODUCTION 


As shoaling waves move from deep to shallow water over a sandy bottom, energy 
dissipation increases in the bottom boundary layer because of the presence of bed 
roughness in the form of wave-induced ripples. Ripple bedforms in coastal shallow water 
play a significant role in the interaction of bottom boundary hydrodynamics and sediment 
transport. When wave energy 1s sufficient, frictional stresses at the bottom cause wave- 
induced ripples to form. The flow over the ripples is dominated by the vortices created at 
the crest on the lee side of the ripples, which also act to enhance wave dissipation. 
Sediment transport is then dominated by avalanches of sand rolling down on the sides of 
the ripples, which is much larger than for a flat bed (Andersen 1999). Wave ripples, 
termed vortex ripples by Bagnold (1946), are triangular in shape with the slope of the 
sides being close to the angle of repose, usually shorter than 60 cm and not higher than 
about 6 cm (Engelund and Fredsoe, 1982). Ripples have complex geometry with 
different equilibrium shapes and sizes. 

Until recently, there have been few observations of bedform generated by waves 
and currents. Some of the first measurements of bedform generation due to waves and 
currents over a sandy bottom were by Inman and Bowen (1963). Tanaka and Shuto 
(1984) performed laboratory experiments on bottom bedforms due to waves and currents 
interacting both in the cross-shore and along-shore directions, and were the first to 
provide an empirical formula for the prediction of ripple geometry due to waves and 


current forcing. 


Most ripple wavelengths fall into three categories of wavelengths (Clifton and 
Dingler, 1984). Orbital scale npples, or “vortex” ripples, have ee (X) that 
scale with wave orbital diameter (2A). Anorbital ripples have wavelengths that depend 
on sediment grain size (d) alone. Suborbital ripples represent a transitionary stage 
between orbital and anorbital ripples and depend on both wave orbital diameter and grain 
size. Important geometric properties to characterize ripples are their height (n), steepness 
(ratio of ripple height to wavelength, 1/4), and cross and along shore symmetry. 
Symmetry in both these directions 1s an important property affecting formation and 
migration rates. A common problem with laboratory data is that the wave periods are 
generally short due to experimental apparatus limitations. Field studies have the 
limitations that the relationship between wave forcing and ripple geometry can be 
complicated due to the range of wave heights, periods, directions, and presence of low 
frequency currents. 

Field experiments have been conducted using a variety of techniques to measure 
bedforms. Manual observations have been made by divers, which greatly limits the 
sample area and sample rate. Boyd et al. (1988) and Amos et al. (1988) used 
photographic measurements to investigate the role of waves and currents in generating 
bedforms. Rotary sidescan sonar was used by Hay and Wilson (1994) and Traykovski et 
al. (1999) to measure the horizontal pattern of bedform evolution. Crawford and Hay 
(1998) developed an underwater video system consisting of a diode laser and a Video 
camera to obtain two-dimensional structures of suspended sediment and bedforms. We 


have adapted a similar approach to analyzing ripple geometry and bedform evolution. 


Measurements described here include waves, ripple formation, and increased 
wave attenuation. This study is a component of the ONR-sponsored os waves 
experiment (SHOWEX) with the overall objective of measuring wave transformation 
mechanisms across the shelf. An underwater video camera in conjunction with a 
structured light source was used to obtain image sequences of the bottom, providing high- 
resolution digital pictures of ripples on the sandy bottom. The ability to resolve 1 mm- 
scale features in the bedforms gives this video imaging system the advantage over 
previous field experiments. The camera and light source were deployed in 13 meter depth 
at the Monterey Inner Shelf Observatory (Stanton 1999). 

Hypotheses for bed migration mechanisms on the inner-shelf (outside the 
surfzone) include wave skewness, mean currents combined with wave stirring, and bound 
long waves. The non-linearity of shoaling waves causes waves to have an onshore 
skewness owing to more peaked crests and flatter troughs. This results in stronger 
onshore velocity and weaker offshore velocity with a net onshore sediment transport. 

Mean currents act to advect sediments that are stirred up and suspended by wave 
action. The direction of transport is in the direction of the current, which could be 
onshore, offshore, or alongshore. 

Bound long waves act to transport sediments offshore. Natural waves are not 
regular, but have a variation in the wave height. A narrow-banded wave train consists of 
wave groups and associated bound long wave (Longuet-Higgins and Stewart 1964). Shi 
and Larsen (1984) recognized that bound long waves under wave groups could give an 
opposing contribution to the transport, because the bound long wave troughs occur under 


high waves with maximum offshore velocity coinciding with high sediment 
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concentrations stirred by the wave group, while the onshore motion occurs under the low 
waves. Deigaard, Jakobsen, and Fredsoe (1999) found the ty of the turbulent 
boundary layer and the coupling between the offshore motion and high waves give an 
offshore mean bed shear stress. From their studies they concluded that bound long waves 
give an offshore contribution to sediment transport. 

The objective of this paper is to analyze observations of bedforms in relation to 
wave and current forcing. The paper is organized into five sections. The second section 
describes the experiment and the equipment that was deployed during September 2 to 
September 24, 1999, calibration and processing data that was acquired. The third section 
describes the results of ripple geometry and ripple migration analyses. A discussion of the 
results is included in the fourth section. The conclusion and recommendations finalize 


this thesis. 


Il. EXPERIMENT 


The MISO instrument frame consists of a five-meter horizontal beam supported 
one meter above the bed by two vertical pipes jetted into the bed as displayed in Figure 1. 
Mounted on the frame is a Broadband Acoustic Doppler Current Profiler (BADCP), 
under-water video camera, paro-scientific pressure sensor, scanning acoustic altimeter, 
Bistatic Coherent Doppler Velocity and Sediment Profiler (BCDVSP) and structured 
laser light. MISO is a component of the Rapid Environmental Assessment Laboratory 
(REAL) and was established July 26, 1999 and continues to be operational. MISO is 
located 400 meters off Del Monte Beach adjacent to the Naval Postgraduate School 
Oceanography laboratory. A local map depicting the location of MISO is shown in 
Figure 2a, and a regional map of California’s central coast is shown in Figure 2b. Waves 
incident to the MISO location in the Southern Monterey Bay tend to be narrow-banded in 
frequency and direction owing to the strong refraction as they pass over the Monterey 
Bay Submarine canyon, and protection by Point Pinos headland. The MISO array is 
located on a sand lens consisting of fine sand with shale outcrops 50 meters on either 
side. The sand at the MISO site was well-sorted quartz with a 0.125 mm median grain 
diameter (d,, ) (Figure 3). 

The data analyzed here is for a 22-day period from September 2, to September 24, 
1999. This period provided comprehensive data from the underwater video camera, 
BADCP, and paro-scientific pressure sensor. The mean water depth was 12.8 m varying 


from 11.7 m to 13.8 m due to diurnal and semi-diurnal tidal forcing. No significant 


storms occurred during this time period. The significant wave height peaked at 1.1 m and 
averaged 0.5 m. | 
A. BROAD BAND ACOUSTIC DOPPLER CURRENT PROFILER (BADCP) 
Three-component velocity profiles with 1-m range resolution were measured 
using the BADCP. The instrument was situated approximately 1 meter off the seafloor 
with the center of the first range bin located 2.1 meters from the transducer face. The 
data were sampled at 0.976 Hz. Spectra were computed using thirty-minute ensembles. 
The spectra were averaged into fifty logarithmic-spaced frequency bins from 0.01 Hz to 
0.35 Hz. Mean-square wave bed velocities were calculated using a linear theory transfer 
function ,( f) to transform measured velocity at 4.1 m, S, ( f), to the bed over a 0.05 
to 0.3 Hz frequency band 


U? - [[H, COT S, COdf (1) 


.05 


where H,(f)- Jost kth 42) k is the wavenumber and h+z is distance above the 


bed. The rms orbital excursion at the bed (A,, ) was computed from the band passed 


brms 


velocity spectrum at the bed using the integrating operator, H ,( f) — Y, , on the 


Am - fa, nr s.Nnar (2) 


To determine the importance of mean currents on ripple formation, thirty-second 


frequency spectrum 


cross-shore and alongshore mean currents were calculated from the BADCP data. It is 


hypothesized that semi-diurnal and diurnal tides, winds and baroclinic tidal effects played 
a role in the forcing of these currents. l 
B. PARO-SCIENTIFIC PRESSURE SENSOR 

A high precission Paro-Scientific pressure sensor was used to estimate surface 
elevation and to determine the mean depth of water. Surface elevation spectra were 


calculated by applying a linear wave transfer function to the pressure spectra, 


cosh(kh) 


S, -[H, C] S, Cf). where pO ERR 


The depth of water was 


calculated from the pressure data averaged every five minutes. The data were acquired at 
2 Hz and with a thirty-minute ensemble. Fifty logarithmically-space frequency bins were 


calculated from 0.001 Hz to 0.2 Hz. Significant wave height was calculated from the 


H, =4 [Sar (3) 


©. VIDEO CAMERA AND LASER LIGHT 


surface elevation spectrum 


Underwater video cameras and structured light sources provide a means to 
observe sediment movement in near shore environments under conditions of good 
visibility. The video/laser system consists of two diode lasers and an imaging camera. 
The diode laser operates in the visible red at 685 nm wavelength with a power of 35 mW. 
The two lasers are housed in a rectangular pressure case 4.6 inches wide by 6 inches long. 
A cylindrical lens in front of each laser spreads the point sources into a 60° beam width, 


4 mm wide sheet of light. One laser sheet of light is directed in the cross-shore direction, 


while the second is directed 45° to the cross-shore direction. The plane of light intersects 
the bed producing a bright stripe in the video image (see Figure 6a). l 

A Cohu 4910, RS 170 monochrome camera is connected to one channel of the 
eight-node power and data controller at the MISO site, which provides serial data 
communication and a video bandwidth channel to shore. This underwater multi-node 
controller provides power to each node and converts the electrical video signal into a 
fiber optic signal for transmission 600 meters to the Marine Operations Building (MOB) 
on the shore. The electrical video signal is converted to a fiber optic signal to transmit the 
long distance while retaining low noise and high bandwidth. At MOB, the signal is 
converted back to a RS 170 signal and digitized by a Haupage frame grabbing board in a 
Pentium III PC running the Linux operating system. The video frames are sampled at 10 
Hz. Each pixel in the images is processed to determine a mean, max, min, standard 
deviation and raw image stored over each five-minute ensemble. Initially the pixel size is 
640 by 480, but is reduced to 320 by 240 in order to conserve storage space and reduce 
noise. 
D. CALIBRATION 

The first step to extract meaningful position data from video images is to calibrate 
the field of view. To calibrate for underwater lens distortion, the camera was suspended 
in a laboratory water tank and exact measurements of distance from camera to a test grid 
were taken. To determine angle versus distance, calibration coefficients were computed 
in the cross-shore axis [cal(x)] and the ripple height axis [cal(z)] as shown in Figure 4. 


The calibrations from pixel space to angle for horizontal angle (Ø) and vertical angle 


(a) is approximately given by 


P = [pix(x) -160]* cal(x) and æ = [pix(z) -120]* cal(z). 
The center pixel in each image was located at x=160 and z= 120 for " 320 by 240 pixel 
sized image. Following the geometry in Figure 5, distance from the camera to the bottom 
is determined by z = d * tan(off + æ), where d is the horizontal distance from the 
structured light to the focal point of the camera and offis the angle between the center 
pixel of the camera and the horizontal frame. For the data considered here, offis equal to 


47.48 degrees and d = 111.76 cm. The cross-shore distance is estimated by 


d 
cc et 
cos(off ) 


Since the instrument frame was tilted 12 degrees in the offshore direction during 


*tan(p). 


these observations, the raw image had to be rotated. Once rotated, the image was 
smoothed by filtering the data using a fourth order low-pass filter with a 2 cm cut off to 
produce a clean line representing relative ripple height in the cross-shore direction, with 
positive values of the cross-section signifying the offshore direction as displayed in 
Figure 6. The low-pass filter had negligible effect on measured slopes. 

During strong wave forcing, turbid conditions existed and video image resolution 
was significantly impaired. However, most of the time visibility was high, and the 
camera provided a high-resolution, two-dimensional time series of the cross-shore bed 


height over a 54 cm cross-shore distance with 1.7 mm/pixel resolution. 


THIS PAGE INTENTIONALLY LEFT BLANK 
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111. RESULTS 


Wave induced ripple formation is related to the ability of waves to move 


sediment, which can be described by the Shields parameter. The Shields parameter is the 
ratio of the force mobilizing the sediments, measured by the bottom stress (7, ), to the 


restoring force of gravity 


Cini (s- Dd © 


where s is the ratio of sediment grain density to water density and is equal to 2.65 for 
quartz sand, g is the force of gravity, and d,, is the mean grain size. Bottom stress can 
be divided into form drag and skin friction, where skin friction is primarily responsible 
for sediment motion and includes the effects of both waves and currents. Assuming a 
classic square law friction factor, and partitioning the stress into wave and current 


contributions, the equation for bed stress is: 


T, =T, tT, == af <ul > +> A «U?; «V!» 6) 


where u,, 1s the rms wave velocity at the bed, U,V are the cross-shore and alongshore 


mean current velocities, and p is water density. The wave friction factor, 
2.5 ds 0.194 e , : 
= le? Gre) “” —5.977] is based on Swart’s (1974) formula using the value 


2.5d,, for the grain roughness of a flat bed (Fredsoe and Deigaard, 1992), where A is 
the orbital excursion radius. The current friction factor is f. =.003 (Amos et al., 1988). 


The initiation of sediment movement occurs when the grain roughness shields parameter 
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is greater than the critical shields parameter 0, = .04 (Grant and Madsen 1976), which 
is consistent with combined mean and oscillatory flows (Amos et al., 1988). 

The evolution of ripple patterns is displayed by examining sequences of the 5- 
minute mean images. Computer animated movies of the mean images provide a detailed 
view of these processes (http://www.oc.nps.navy.mil/~stanton/ miso/dixon). The 22-day 
period of observation from yearday 245 through yearday 267 coincided with two wave 
events, two to three days in duration (Figure 7a). A wave event also preceded yearday 
245, which is included in the wave climatology. The wave-induced cross-shore bed 
velocity (Figure 7b) reflects the changes in wave height. The mean currents (Figure 7c) 
were weak with alongshore currents (<0.1m/s) stronger than the cross-shore currents 
(<0.05m/s). Diurnal variation in these currents (either due to tide or diurnal sea breeze) is 
evident. The combination of larger wave events succeeded by relatively small wave 
events created highly variable bedform evolution. 

Ripple wavelengths and migration are summarized in Figure 8a, which shows 
ripple height (grayscale intensity) versus cross-shore distance and time. The dotted lines 
are the crests of two ripples determined by a change in slope from positive to negative. 
Migration rates are determined from the slope of the crest line. The table describes nipple 


crest migration rates. 


246-248 248.5-249.4 | 254-257 259.5-262 | 264-265.3 
Ehe 





Crest#2(cm/day) | 10.35 -6.23/7.07 LAT 


Significant migration occurs between 246-248, 248.5-249.4 and 259.5-262. On yearday 


248.5, a bifurcation occurs as crest #2 separates, leaving one crest that continues to 


migrate offshore, while a second crest migrates onshore. On yearday 265, it is interesting 
to note that while forcing is low (0 <0,,, ), crest +1 migrates offshore fairly swiftly but 


crest #2 does not move. It is not understood what causes these differences in migration 
directions, however it is hypothesized that the crests are associated with a more three 
dimensional bedform. 

The Shields parameter was calculated using the measured BADCP velocities 
partitioned into wave mean square velocities (Equation 1) and mean current velocities. 
The wave velocities were approximately an order of magnitude greater than the mean 
current velocities. In addition, the wave friction factor is an order magnitude greater than 
the current friction factor. Hence, the contribution to the Shields parameter by the mean 
current was negligible (2 percent on average). 

Marginally active ripple patterns and evolution were seen throughout the period 


of observation despite the fact that the grain roughness Shields parameter (Equation 4) 


often dips below the critical Shields parameter 8 


e = -04 (Figure 8b). The video imaging 
system did not produce bedform images during daylight hours because the solar 
illumination usually swamped the structured laser light sheet. Therefore, each day is 
missing approximately 7 hours of images. Images are also missing during a large wave 
event between yeardays 249.5 and 252, because turbidity is strong during high sediment 
suspension. 


Ripple length ZA , is the distance from crest to crest and ripple height, 77, is the 


difference between the elevations of the crest and succeeding trough. A maximum angle 
13 


($) of the ripple is calculated as the maximum angle occurring between a crest and 


succeeding trough, while the slope angle 1s the average angle between the crest and 


succeeding trough, VA . Maximum angle ($) was determined by computing the difference 


in tan” Cos ) with Ax = 0.8 cm being the image resolution. The slope angles were 


found not to be sensitive to the choice of Ax values, which was varied between 0.8 and 
2.4 cm. 

Images of the evolving ripple geometries are shown in Figure 9, where 3-4 day 
time series of data are plotted against relative ripple height and cross-shore distance. 
Observations begin with relic npples during yearday 245, which are left over from the 


large wave days between yeardays 242-244. Relic ripples have smaller slope angles and 
more rounded crest than active vortex nipples. By yearday 246, @ exceeds @.,, and 
wave ripples begin to form. These are pointed, active vortex ripples, which shed vortices 


regularly. The Shields parameter remains above the 8 


er threshold between yeardays 246 
and 249, during which time ripple migration is weakly offshore; however, due to the lack 
of bed profiles during these days, the ripple migration is difficult to resolve. The ripples 
are moderately asymmetric containing a steeper slope on the leeward (off-shore) side. 
During yearday 249, strong wave forcing caused strong sediment suspension, which 
blocked the vision of the video camera so that bed profiles were not recovered between 
yeardays 249.5 and 252. During this two and a half-day period, the Shields parameter 


reached a maximum over the twenty-two day time period. It cannot be determined if the 


ripples migrated onshore or offshore during this period of time 


14 


On yearday 252, the waves calmed with@ below 6 _, until yearday 254 leaving 


crit 
relic ripples, which significantly increase in size as only two crests dominate over the 57 


cm cross section. The ripple length of 30 cm between these successive crests is the 


longest during the entire period. Between yearday 252-254, the ripples are inactive 
representing relic ripples. From yearday 254-256, 0 is above the threshold value and 


bedforms begin to move slowly, migrating at a rate of 4-5 cm/day. This well organized 


long crested wave ripple is sustained as the bedform migrates seaward. 


Ripple migration ceases on yearday 256.5 as @ is below 0 


crit 


and relic ripples 
prevail. The long crested wave ripple remains intact until yearday 258.5 when the second 


largest wave day begins with @ well above 8 


«i» aNd the ripples are active. Vortex ripples 
occur during the larger wave days, and their geometry is highly variable. Ripple 
migration begins at the same time that the large swell arrives (yearday 259). Migration 
rates are equal to that of the previous period at 5 cm/day. The long crested ripples that 
dominated earlier give way to shorter bifurcated three-dimensional ripples during 
yeardays 261 —263 (Figures 10a and 10b). Two ripple crests dominated the bed between 
252-259. The wave event (yearday 259) stirs the bed producing the bifurcation and three- 
dimensional ripples that are seen on yearday 261. Between yearday 261 and 265, wave 
energy is weak, and ripples are less active but continue to migrate a few cm/day. Bed 
movement slowly decreases and ripple migration comes to a halt on yearday 264. There 


is a small burst of wave energy on yearday 265, and ripples are highly active. From 


yeardays 266-267, wave energy is low and bedform variability is limited. The three- 
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dimensional bifurcated rıpple gives way to symmetric gently sloping ripples on yearday 
267. l 

In summary, the observation period is marked by three wave events. The first 
occurs prior to the start of ripple observation, the second wave event (249-252) is too 
strong to obtain bedform images, and bedform geometry is highly variable during the 
third wave event. It is during these larger wave events that ripple migration rates are the 
highest while ripples are strongly asymmetric with steeper slopes on the leeward side of 


the crest. Immediately following the larger wave events, relic ripples prevail as@ is 


below 0... Once wave energy increases slightly, ripples begin to migrate, although their 


crit ° 
geometry does not significantly change. The migrating vortex ripples yield to relic 


ripples as wave energy decreases and the pattern starts over. 
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IV. DISCUSSION 


Ripple steepness is limited by the angle of repose of the wet sand, where ø is the 
local maximum angle between crest and succeeding trough. Assuming the ripple is a 


triangle, Nielsen (1981) suggests that the maximum steepness a bedform can obtain 
(v, ) =0.5tand. Ifthe trough had the shape of a parabola, the limit would be 


(7/A) =0.25tan @ as displayed in Figure 11. Most of the bed slope data from the 
observations lies within the interval 0.25 tang <(7/A)<0.S5tang (Figure 12). Nielsen 
(1979) found that ripples maintain their maximum steepness when @ < 0.2 but, when 


0>0.2, the vortex can no longer maintain the maximum steepness against the erosion of 


the crest. The result is that 77/ A decreases as water particle velocity increases, and at 
stronger forcing ripples disappear completely. During the observation period,@ peaked 
at 0.12 and we do not see 7/A decreasing as a function of water particle velocity (Figure 
indicative 


13a). Ripple steepness varied between 0.04 and 0.26 at low values of 0 «0 


of relic ripples. 
The ripple length has been suggested to be related to the orbital diameter (2A) of 


the wave-induced fluid motion near the bed, 4 z 1.34 (Nielsen 1992). However, we find 
significant variation of over the range 0 « 0.05 , with ^ decreasing to an 
approximate value of 1.3 (Figure 13b). One has to be careful in the interpretation of the 


apparent decrease in Y, with increasing mobility number as 0 « u* = A 
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Therefore, both the ordinate and abscissa contain the same variable. Similar plots have 
been widely reported in the literature (Nielsen, 1992). For the larger value of 6 , un 


approaches 1.3 independent of @ in the range of 9 >0.11 up to the maximum value of 


0.21. 


Clifton and Dingler (1984) summarized several data sets from laboratory and field 


studies and suggested that that ripple wavelength scales linearly with orbital diameter for 


shorter period waves until orbital diameter/grain size Q4, ) reaches a value of 2000. 
50 


These are described as orbital or vortex ripples (Bagnold, 1946). In the range of 


2000« ZM <5000, for longer period waves with increased orbital excursion, the ripples 
50 


go through a transitional stage, where ripples are much shorter and are referred to as 


suborbital ripples. For values of a >5000, when orbital diameter is larger, ripple 
50 


wavelength scales directly with grain size. The bedform environment was in a state of 
flux during the observation period. Orbital and suborbital ripples were the dominant 
ripple pattern (Figure 14). However, when wave forcing was small, relic ripples were 


present. The period starts with remnant or relic ripples, which are left over from the 


larger waves on yearday 243-244. These appear to be suborbital relic ripples as — 
50 


fell between 2000 and 5000 on yearday 244. Ripples tend to orbital ripples on yearday 
246 as ripple length scales approach orbital excursion. The bed is in constant motion 
with ripple steepness peaking on yearday 246 and the steepness of these vortex ripples 


approach their angle of repose. When the largest waves of the observation period arrive 
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on yearday 249.4, no images are retrieved, but it can be hypothesized that these are 
highly active suborbital ripples. When the waves subsided, suborbital rei ripples were 
left. The strong wave forcing changed the bedforms to two ripple crest spaced 30 cm 
apart. When the next large wave event occurs on yearday 259, ripple length scales with 
orbital excursion and orbital ripples are formed. Once the forcing ceases, orbital relic 
ripples remain and continue through to the end of the observation period. 

Ripples become asymmetric when suspended sediments are lifted and carried 
from the upstream side trough and settle on downstream side crest. In the following 
discussion, upstream is the side of the ripple where the sediment is moving from and 
downstream is the location where the sediment is moving. The building up of the 
downstream side crest creates a steeper slope on the downstream side and a shallower 
slope on the upstream side. As this process continues, the downstream side becomes 
steeper and the upstream side becomes shallower. Active asymmetric ripples migrate in 
the direction of the steeper slope, while active symmetric ripples may not change their 
position. For example, during yearday 254, the leeward slope is twice as steep as the 
windward slope indicative that the ripples are migrating in the offshore direction (see 
histograms of slope angle for selected ripple profiles in Figure 15). 

Mean slope angles were calculated by subtracting the mean slopes of the offshore 
and onshore sides of the ripple to determine ripple symmetry. Large positive values 
represent asymmetric ripples with a steeper slope on the offshore side, indicative of 
ripples migrating offshore (negative slopes indicate onshore migration) while small 
values correspond to symmetric ripples (Figure 16). Asymmetric ripples were formed 


following the wave events and slowly eroded with weaker forcing. The forcing by larger 
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waves on yeardays 243-244 appears to have resulted ın the observed asymmetric ripples 
on yeardays 245-246. These bedforms became symmetric as the — forcing subsided. 
The larger waves starting on yearday 248 had large onshore, negative skewness (see 
Figure 17) coinciding with onshore migration (see Figure 8a) and ripple asymmetry 
indicative of onshore transport. From yeardays 250-253, @ decreased and skewness 
changed from onshore to offshore. When the images resumed, ripples migrated onshore 
(yeardays 253-257) and were asymmetric. As wave forcing decreased, migration rates 
slowed and ripples became symmetric. During the last wave event (yeardays 259-262) 
migration rates were strong. Ripples became asymmetric for a short time on yearday 263 
and tended toward symmetry as the period of observation concluded. 

Mean and skewness velocities were calculated to assess if they could explain 
offshore ripple migration. Skewness of low-passed and high-passed filtered velocities 
(0.05 Hz cut-off), most of the time showed no preferred direction (Figure 17). The 
exceptional times when onshore, negative skewness exceeded 0.1 occurred on yeardays 
250 and 265 when migration was onshore. Cross-shore mean current were weak and 
fluctuated between the offshore and onshore direction, again suggesting no preferred 
direction (see Figure 7c). 

Bound long waves give an offshore contribution to sediment transport and can 
result in offshore npple migration as schematically shown in Figure 18. The presence of 
bound long waves was measured by correlating the infragravity wave band (obtained by 
low passing the cross-shore velocity below 0.05 Hz) with the envelope of the sea-swell 
band of frequencies (Hanes 1991) as shown in Figure 19a. The envelope of the sea-swell 
band is obtained by low pass filtering (below 0.05 Hz) the u* time series. Positive 
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correlation suggests bound long waves, (coordinate system is u positive offshore). The 
correlations are significantly different than zero at the 95% confidence level (assuming 
Guassian velocities) for -0.0167< r <0.0167. The r values are persistently positive and 


mostly significant indicating bound long waves were present 


A measure of sediment transport advected by bound long waves is given by 
u,” (u?) p Where sediment transport is parameterized as proportional to (u?) ¡p (Figure 
19b). Transport was mostly offshore in agreement with ripple migration direction. 


Four offshore ripple migration events occurred between 246-248, 254-257, 259.5- 


262, and 264-265.3 as shown in Figure 9. The first event yielded the swiftest migration 


rates. Orbital forcing was strong enough for@ 


crit 


threshold to be surpassed. Two ripples 


dominated at the start of yearday 246 and the ripple #2 migrated at 10.35 cm/day while 
ripple #1 migrated at 6.84 cm/day. During this time frame, waves were groupy and bound 
long waves were present. The migration rate was low during the second event as O, , was 
exceeded only during yearday 255 by 0.02. These ripples migrated offshore at 1.5 
cm/day. Bound long waves were present and ripples slowly migrated offshore. Yeardays 
259.5 and 262 provide good information because although the forcing was strong, video 
images were obtained. The presence of bound long waves caused the ripples to migrate 
offshore at 6.3 cm/day. 


There 1s a question of what causes the final event to migrate offshore at 5.0 


cm/sec when @ < @ 


si- This brings into question the validity of 9_,=0.04 or how this 


ril 


value is applied. By viewing the video showing ripple forms migrating along with 


calculated Shields parameter, it is qualitatively observed that the bed only migrated when 
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0 >06 with exception of yearday 264. The basis for@.,, =0.04 are laboratory 


experiments using monochromatic waves with constant amplitude (Grant and Madsen, 
1976). The wave velocity is a bulk rms value for a 30 minute ensemble and does not 
accurately capture the large range of forcing between the peak and minima within a wave 
group. For very narrow banded, groupy waves in nature, the maximum wave velocity 


could be more two-five times the monochromatic value for the same Shields Parameter. 


Therefore, it is hypothesized 6_,,, may be less than 0.4 for waves in nature. 
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V. CONCLUSION 
The analysis of bedform evolution using an underwater video camera in 
conjunction with a laser light has been demonstrated. High-resolution images were 


obtained and analyzed in relation to wave and current forcing. Ripple geometry and 


migration rates were obtained. The following conclusions are made. 


Ripple steepness varied for values of @ , and no relationship was found. Relative 


ripple length Y, ) decreases as a function of the Shields parameter and was 


approximately constant at 1.3 in the range of 0.07 «0 « 0.21. Ripples were steeper on 
the downstream or offshore side. Asymmetric ripples were left after large wave events 
and the tended toward symmetry with weaker forcing. Orbital ripples were the dominant 
bedform found when the wave and current forcing was strong. However, when the waves 
relaxed relic ripples prevailed. Forcing was never strong enough such that orbital 
excursion could reach a value to create anorbital ripples, however, suborbital ripples were 
formed during the strongest forcing event. 

Migration of ripples may be the dominant mode of cross-shore sediment transport 
at the MISO site. Measured migration rates occurred in four events (246-248, 254-257, 


259.5-262, and 264-265.3), and all ripples actively migrated offshore. During these 
events, O exceeded the critical shields parameter on all days except during yearday 264- 


265.3. Skewness of orbital velocity (both high and low passed) mostly showed no 
preferred direction to explain offshore ripple migration. Significant positive correlation 


was obtained between wave groupiness of the narrow banded swell, (measured as the 
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envelope of high passed velocity signal), and bound long waves, which may account for 
the offshore migration during the first three events. 
The use of video imaging to accurately map the bedform evolution was 


successful, however when large waves caused strong sediment mixing, no video images 


were obtained. Hence ripple information was limited to low wave events. 
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STRUCTURED LIGHT 


mel! 





M 





JETTED PIPE 


Figure 1. Instrument frame (viewed from sea to shore) at MISO array. ADCP was used to 
measure current velocities. The PARO was used to calculate wave height and water 
depth. The camera and structured light provided high resolution video images of 
bedforms. 
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Figure 2a. Location of MISO instrument array. The star represents the location of the 
MISO array. Figure 2b. Regional map of central coast of California showing that the 
Point Pinos headland protects the MISO location and causes waves to be narrow-banded 
in frequency and direction owing to strong refraction. 
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Sand Sediment Distribution at MISO Site - 
0.4 


0.3 


Sand Fraction 
C 
ho 


Sieve Size (mm) 


Figure 3. Sand Grain Size Distribution with median grain size (d,, )=.125 mm 
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Actual Angle (Radians) vs Horizontal Distance (Pixels) 


equation of line is y= 0.00051656 x + -0.00025728 
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Figure 4a. Horizontal calibration Figure 4b. Vertical calibration. The slope of the line is 
the calibration coefficient. 
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Figure 5. Camera calibration showing horizontal and vertical cross-sections. Alpha and 
B are angles in the vertical (z) and cross-shore (y) directions. 
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Figure 6a. Raw image taken from video camera. Figure 6b. Conversion of video image 
in pixel space to raw image in real space (cm). Figure 6c. Image was rotated because the 
instrument frame was tilted 12 degrees offshore. Figure 6d. The image was then 
smoothed by applying a low pass filter with a 3cm scale. The filter did not affect the 
image resolution which was 1.7mm/pixel. 
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Figure 7a. Significant Wave Height ( H, ) for yearday 243-267 of 30-minute averages. 
Figure 7b. Rms orbital velocity at the bed (u, ) of 30-minute averages. Figure 7c. 
Mean currents (cross-shore and alongshore of 1-hour averages). 


3] 


cross -shore distance (cm) 


Shields Parameter 


Figure 8a. Time series of Relative npple height displayed in grayscale (white is large, 
dark is small) vs. cross-shore distance (positive values are offshore) and time. Ripple 
Crest #1 1s the lower npple crest and Ripple Crest #2 is the upper ripple crest. Greatest 
migration rates occur between 244-248. Figure 8b. Grain roughness Shields parameter 
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calculated from a rms orbital velocity averaged every 30 minutes. 
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Figure 9. Time evolution of ripple geometries displayed in 3-4 day periods of relative 
ripple height (grayscale) vs. cross-shore distance (positive is offshore). 
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Figure 10a. Raw image of long crested ripples on yearday 253 (9/10/99) evolve into 
shorter bifurcated ripples (Figure 10b) on yearday 261 (9/18/99). The long crested two 


are left after the strong forcing on yearday 249-252. These give way after the high waves 
on yearday 259-261 to a short crested bifurcation. 
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Figure 11. With maximum slope on the npple profile is equal to tan@, the height to 
length ratio must for geometric reasons be 0.5 tan @ for a triangle npple and 0.25 tang 
for a parabolic ripple (after Nielsen 1992). Í 
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Figure 12. Time series of ripple steepness/tan 9. Bed slope geometry data generally falls 
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Figure 13a. Ripple steepness is highly variable. Figure 13b. Relative npple length 
parameterized by orbital excursion (plotted only when 0 > 0... ) decreases as forcing and 


crit 


orbital excursion increases. The dashed line corresponds to 4 z 1.3A. 


an 


Ripple Geometry 


Anorbital 


Suborbital 





245 250 255 260 265 


gI 
oO 


I 
| 
I f 
I 


A 
o 


RemnantOrbitalNo Data ; Remnant orbital Remnant 
| Suborbital | ^. Orbital 


- 


e 


N Q) 
O oO 


Ripple Length & Orbital Excursion 
E 





245 250 259 260 265 
yearday 


Figure 14a. Ripple Geometry (Clifton and Dingler, 1976) classified into orbital, 
ed ) < 2000, suborbital, 2000 < M < 5000, and, anorbital, Y >5000. Figure 
50 50 50 


14b. Time series of A and A. Only values of A have been plotted when @ > 8m- 
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Figure 15. Left panel show evolution of selected images of cross-shore distance (positive 
is offshore) vs relative npple height. Right panel shows histogram of slope angle with 
mean slope angle displayed above histogram. Mean slope angles were calculated by 
averaging the slope angle between crest and succeeding trough. The negative slopes 
angle are the mean slopes on the off shore side, while the positive slope angle are the 
slope angles on the onshore side. 
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Figure 16a. Ripple Asymmetry: Large positive values show steeper sloping rıpples on the 
off-shore side indicative of offshore migration while small values show symmetric 


ripples. Negative values show steeper sloping mpples on the on-shore side. Asymmetric 


ripples were left after strong wave forcing (high value of Shields Parameter, Figure 16b) 


and tended toward symmetry as the forcing subsided. 
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Figure 17. Skewness of cross-shore orbital velocity from surface gravity waves. A high 
pass filter at 0.05 Hz was applied to the cross-shore velocity time series and skewness 
values were calculated over one-hour samples. 
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BOUND LONG WAVE THEORY 
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Figure 18. Change in momentum (high waves) is balanced by set down of forced long 
wave and small waves is balanced by a relative set-up. The forced long wave gives 
reverse flow under high wave and forward flow under small waves. Because more 
sediment is suspended during the larger waves. Net transport is revered or offshore. 
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Figure 19a. Correlation Coefficient determined by correlating the infragravity wave band 
(obtained by low pass filtering the cross-shore velocity at 0.05 Hz) with the envelope of 
the sea-swell band of frequencies. The dashed lines represent the significant correlation. 
Figure 19b. Correlation as calculated in (a), but is not normalized. Correlation is 


determined by u, * (u^),,. 
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